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Abstract 

Niemann-Pick Type C (NPC) disease is an autosomal recessive neurodegenerative disorder caused in most cases by 
mutations in the NPC1 gene. NPC1 -deficiency is characterized by late endosomal accumulation of cholesterol, impaired 
cholesterol homeostasis, and a broad range of other cellular abnormalities. Although neuronal abnormalities and glial 
activation are observed in nearly all areas of the brain, the most severe consequence of NPC1 -deficiency is a near complete 
loss of Purkinje neurons in the cerebellum. The link between cholesterol trafficking and NPC pathogenesis is not yet clear; 
however, increased oxidative stress in symptomatic NPC disease, increases in mitochondrial cholesterol, and alterations in 
autophagy/mitophagy suggest that mitochondria play a role in NPC disease pathology. Alterations in mitochondrial 
function affect energy and neurotransmitter metabolism, and are particularly harmful to the central nervous system. To 
investigate early metabolic alterations that could affect NPC disease progression, we performed metabolomics analyses of 
different brain regions from age-matched wildtype and Npcf'' mice at pre-symptomatic, early symptomatic and late stage 
disease by ^-NMR spectroscopy. Metabolic profiling revealed markedly increased lactate and decreased acetate/acetyl-CoA 
levels in NpcV 1 ' cerebellum and cerebral cortex at all ages. Protein and gene expression analyses indicated a pre- 
symptomatic deficiency in the oxidative decarboxylation of pyruvate to acetyl-CoA, and an upregulation of glycolytic gene 
expression at the early symptomatic stage. We also observed a pre-symptomatic increase in several indicators of oxidative 
stress and antioxidant response systems in NpcT 1 ' cerebellum. Our findings suggest that energy metabolism and oxidative 
stress may present additional therapeutic targets in NPC disease, especially if intervention can be started at an early stage of 
the disease. 
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Introduction 

Niemann-Pick Type C (NPC) disease is a fatal, autosomal 
recessive neurodegenerative disease caused in 95% of cases by 
mutations in the NPC1 gene and in the remaining cases by 
mutations in JVPC2. NPC1 is a late endo/lysosomal protein that 
interacts with NPC 2 in the lysosome lumen to mediate cholesterol 
egress from endosomes to plasma membrane and endoplasmic 
reticulum [1-3]. The loss of either NPC1 or NPC 2 leads to the 
endo/lysosomal accumulation of unesterified cholesterol, an 
impaired cholesterol homeostatic response and a wide range of 
other cellular abnormalities, such as decreased oxysterol produc- 
tion, alterations in calcium homeostasis and autophagy, increased 
oxidative stress and inflammatory responses [4—9]. The most 
severe consequences of NPC 1 -deficiency are a progressive loss of 
Purkinje neurons in the cerebellum and widespread neuronal 
abnormalities. Other common characteristics of NPC 1 -deficient 
brain include altered synaptic transmission, glial activation, 



neuroinflammation, and dysmyelination [10-13]. In the periph- 
ery, increased total cholesterol levels cause hepatosplenomegaly 
and liver dysfunction [14]. Neuronal cholesterol accumulation 
begins at early, pre-symptomatic stages of the disease, and 
mobilization of endosomal cholesterol in vitro or in vivo can prevent 
many of the cellular abnormalities and symptoms of NPC disease 
[15-19]. Although the defect in cholesterol trafficking causes 
many, if not all, cellular abnormalities observed in NPC disease, 
relatively little is known about the cellular mechanisms leading to 
neuronal dysfunction and neurodegeneration. 

Endosomal cholesterol is transported to mitochondria even in 
the absence of functional NPC 1 , and the endosomal cholesterol 
accumulation can lead to increased cholesterol levels in mito- 
chondria [20,21]. Increased levels of mitochondrial cholesterol 
may lead to lower ATP production and a greater sensitivity of 
these mitochondria to oxidative stress [22,23]. In human 
embryonic stem cell-derived NPC 1 -deficient neurons, defective 
mitophagy led to accumulation of depolarized, fragmented 
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mitochondria [24] . Several reports of increased oxidative stress in 
NPC 1 -deficiency further implicate mitochondria in NPC disease 
pathology. For instance, NPC patients have increased serum levels 
of oxidative stress markers [8,25,26], and oxidative tissue damage 
was observed in the liver and cerebellum of late symptomatic 
NPC 1 -deficient mice together with a gene expression profile 
indicative of increased oxidative stress [9,27]. 

The brain is particularly sensitive to alterations in mitochondrial 
integrity or function due to its high energy requirement and 
reliance on oxidative metabolism [28]. Most neurodegenerative 
diseases are associated with mitochondrial dysfunction and 
primary mitochondrial defects commonly manifest in the brain 
[29]. The compartmentation of pathways between neurons and 
astrocytes and the close connection between energy metabolism 
and neurotransmitter balance further increase the complexity of 
brain metabolism [30-32]. 

Here we have used an unbiased metabolomics approach to 
investigate changes in brain energy metabolism beginning at pre- 
symptomatic stages of NPC disease. Metabolites were measured by 
'H-NMR spectroscopy in the cerebella, cerebral cortices and 
hippocampi of BALB/cNctr-Npcl mlN /J mice, a well-character- 
ized murine model of NPC disease with a null mutation in Npcl. 
Homozygous NpcT 1 ' mice are asymptomatic at birth, but develop 
progressive ataxia, tremors and hindlimb dystonia around 5-6 
weeks of age, and have an average lifespan of approximately 1 2 
weeks of age [33]. The pathological changes in the brain of these 
mice at different stages of the disease have been described in some 
detail [10,33-36]. Protein and mRNA expression analyses of key 
metabolic proteins gave further insight into mitochondrial changes 
and into the roles of neurons and astrocytes in NPC 1 -deficient 
cerebellum. Our results revealed significant pre-symptomatic 
alterations in glucose and pyruvate metabolism, as well as an 
activation of stress-related pathways in NPC 1 -deficient brain. 
Further investigations will be required to determine whether the 
identified pathways could present additional targets for therapeutic 
intervention to alleviate the symptoms of NPC disease. 

Materials and Methods 

Materials 

Neurobasal medium, Minimum essential medium (MEM), B27 
and fetal bovine serum were purchased from Life Technologies 
(Burlington, ON). Poly-D-lysine was obtained from Peptides 
International (Louisville, KY). [ 3 H] -deoxyglucose (2-[l,2- 3 H(N)]; 
1 mCi/ml; 25-50 Ci/mmol) was obtained from Perkin-Elmer 
(Waltham, MA). Other chemicals and buffers were obtained from 
Sigma (Oakville, ON) or Thermo Fisher Scientific (Ottawa, ON). 

Animals 

NPC 1 -deficient mice were derived from an in-house breeding 
colony of BALB/cNctr-Npcl mil VJ mice heterozygous for a null 
mutation in Npcl (originally obtained from The Jackson Labora- 
tory, Bar Harbour, ME, strain 003092). Only wildtype mice (WT) 
and mice homozygous for the mutation {Npcl~ / ~) were used in this 
study. All animals were genotyped by PCR with primers as 
described [34]. Mice were fed a normal chow diet (Prolab RMH 
3000 5P00*, LabDiet) ad libitum and not fasted prior to 
termination. To collect brain tissue for analysis by 'H-NMR or 
quantitative PCR (qPCR), animals were deeply anesthetized with 
halothane (Sigma, Oakville, ON) and terminated by decapitation 
at postnatal day P21/22, P35/36, or P49/50. Brain tissue was 
removed rapidly, and the cerebellum, cerebral cortices and 
hippocampi were dissected and immediately snap-frozen in liquid 
nitrogen. 



Ethics statement 

All procedures were approved by the animal ethics committee 
of Dalhousie University (protocol numbers #10-102 and #10- 
030) based on the standards established by the Canadian Council 
of Animal Care. 

1 H-NMR spectroscopy 

Metabolites were extracted for 1 H-NMR spectroscopy essen- 
tially as described [37]. Briefly, tissues were mixed with methanol 
to a final methanol/water ratio of 3.2:1 and homogenized via 
ultra-sonication in a FastPrep®-24 instrument (MP Biomedicals 
Inc., Solon, OH). Homogenates were extracted with chloroform 
and water in a final chloroform/ methanol/ water ratio of 1 : 1 :0.8 1 . 
Aqueous extracts were dried under nitrogen, re-dissolved in 
700 uL of phosphate buffer (1 mM NaH 2 P0 4 , pH 7.2) with 
172.2 mg/mL of sodium 3-trimethylsilyl-2,2,3,3-d 4 -propionate 
(TMSP, Cambridge Isotope Laboratories, Andover, MA) in 
D z O as internal reference standard, and stored in liquid nitrogen 
until analysis. lD-'H-NMR spectra were acquired using a 5 mm 
TCI CryoProbe™ (Bruker Biospin) on a Bruker Avance III 
spectrometer (Bruker Biospin) operating at 700 MHz proton 
resonance frequency [38] . A detailed description of the acquisition 
parameters and conditions is attached as Text SI. 

Gene expression analysis by quantitative PCR 

Total RNA was prepared from snap-frozen brain tissue using a 
commercially available kit (Aurum Total RNA kit, BioRad, 
Mississauga, ON). One (J.g of RNA was used to prepare cDNA 
using iScript reverse transcriptase (BioRad) according to the 
manufacturer's instructions. Quantitative PCR (qPCR) was 
performed on 0.66 ng of reverse transcribed RNA using the 
iSYBR Green Mastermix (BioRad) and primers as listed in Table 
SI. Specificity of the amplification reaction was tested by melt 
curve analysis, agarose gel electrophoresis, and sequencing of the 
amplicon. Data were calculated by the Pfaffl method with 
cyclophilin (Ppia) as a housekeeping gene and standardized to 
the wildtype samples of the same age and the same run of the 
qPCR analysis [39]. Two additional housekeeping genes, fi-Actin 
and Rpll3a were also tested. Calculation of Rpll3a expression 
using either Ppia or f}-Actin as a housekeeping gene yielded marked 
differences between wildtype and NpcT*' mice at 3 weeks of age 
(Table S2), making Rpll3a unsuitable as a housekeeping gene. Ppia 
and fi-Actin were not differentially expressed between wildtype and 
NpcT 1 ' mice when calculated relative to each other (Table S2), and 
gave similar results when used as housekeeping genes for 
expression analyses of target genes. All qPCR data shown in the 
main body of the manuscript were standardized to Ppia. Data were 
derived from 8 to 15 mice of each genotype and each age. 

Immunoblotting 

Snap-frozen tissue was homogenized in ice-cold HEPES buffer 
(10 mM HEPES, 1 mM EDTA, 1 mM EOT A, 1% Triton X-100, 
0.5% Nonidet 40), with freshly added protease inhibitors (5 ug/ml 
leupeptin, 5 Ug/ml aprotinin, 50 |J,M PMSF and 1 uM pepstatin) 
and phosphatase inhibitors (2 mM ortho-vanadate and 1 mM 
sodium fluoride) using a motor-driven pestle fitting into a 
microcentrifuge tube in a metal cooling block, followed by slow 
trituration through a syringe fitted with a 26 gauge needle. Protein 
content was determined by a bicinchoninic acid-based photomet- 
ric assay (Thermo Fisher Scientific). Brain homogenates were 
separated by reducing SDS-PAGE and transferred to polyvinyli- 
dene fluoride membranes. The membranes were blocked in 5% 
skim milk powder in Tris-buffered saline with 5% Tween (TTBS, 
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pH 7.4) and incubated with primary antibodies as indicated. 
When anti-phospho antibodies were used, membranes were 
blocked with TTBS containing 2% polyvinylpyrrolidone and the 
phosphatase inhibitors sodium orthovanadate and sodium fluo- 
ride. Secondary horseradish peroxidase-conjugated donkey anti- 
mouse, anti-rabbit, anti-chicken, or anti-goat antibodies (Jackson 
Immunoresearch, West Grove, PA) were diluted 1:10,000 in 
TTBS, and detected by enhanced chemiluminescence. Actin was 
used as a loading control. Primary antibodies are listed in Table 
S3. Immunoblots were performed with samples from at least 6 
mice of each genotype and age. 

Determination of mitochondrial DNA 

DNA was isolated from 5-week old wildtype and NpcT 1 ' 
cerebellum by phenol-chloroform-isoamyl extraction as described 
[40]. The cytochrome c oxidase subunit I (Mtcol) gene of 
mitochondrial DNA and the JVdufil gene of nuclear DNA were 
amplified by qPCR of 0.1 ng DNA using SYBR Green and 
primers as described [41] and listed in Table SI. Specificity of the 
amplification reaction was tested by melt curve analysis, agarose 
gel electrophoresis, and sequencing of the amplicon. Data were 
calculated by the Pfaffl method using Ndufvl as standard and Mtcol 
as target gene [39]. Data were derived from 10 mice of each 
genotype. 

ATP generation by isolated mitochondria 

Cerebella, cerebral cortices and hippocampi were dissected 
from 3- or 7-week old wildtype and NpcT mice and homogenized 
together in cold mitochondria isolation buffer (220 mM mannitol, 
7 mM sucrose, 20 mM Hepes, pH 7.2, 1 mM EGTA, 0.1% 
bovine serum albumin (Fraction V, fatty acid free) and protease 
inhibitors. Mitochondria were isolated by differential centrifuga- 
tion. ATP production was detected by the luciferase-luciferin 
method as described [42]. Briefly, mitochondria (100 or 50 (ig) 
were incubated with energy substrates (1 mM pyruvate/ 1 mM 
malate or 5 mM succinate), ATP detection buffer (40 \lM 
luciferin, 1/3000 luciferase (Biotium)) and ADP (300 uM) in Tris 
buffer (150 mM KC1, 25 mM Tris pH 7.4, 2 mM EDTA, 0.1% 
fatty acid free BSA, 10 mM KH 2 P0 4 , and 0.1 mM MgCl 2 ) at 
37°C. Luminescence was measured with a FLUOStar Optima 
plate reader (BMG Labtechnologies) for 20 readings at 1 s 
integration time. The rate of ATP generation was calculated as 
the luminescence increase during the measurement period of 20 s 
as luminescence (a.u.) per second per cell protein, and expressed as 
a percent of the average of wildtype samples of the same 
experiment. Addition of NaNs (10 mM), rotenone (2 |J,g/ml), or 
antimycin A (50 ug/ ml) reduced the luminescence increase during 
the measuring period to undetectable levels even in the presence of 
energy substrates (not shown). To verify mitochondrial polariza- 
tion, an aliquot of each mitochondrial preparation was incubated 
with 2 uM nonyl acridine orange (Biotium, Hayward, CA) for 
30 min at 37°C, washed, and analyzed for fluorescence at 488/ 
520 nm excitation/emission. 

Neuronal glucose uptake 

Primary cortical neurons were prepared from embryonic El 7 
wildtype and NpcT 1 ' mice essentially as described [43]. Neurons 
were grown in serum-free Neurobasal/MEM medium (1:1 v/v), 
with 2% B27 supplement, 0.25 |0,M glutamate, and antibiotics for 
10 days. Five days after plating, neurons were treated with 2.5 uM 
cytosine arabinofuranoside to prevent proliferation of glia cells. 
For measurement of glucose uptake, neurons were incubated with 
osmolarity-adjusted HEPES-buffered saline (10 mM HEPES, 
124 mM NaCl, 3 mM KC1, 2 mM CaCl 2 , 1 mM MgCl 2 ; 



osmolarity 270 to 290 mOsm), containing 2 mM glucose and 
2 |iCi/ml [ 3 FTJ-deoxyglucose for 30 min. Cell-associated radioac- 
tivity was determined by scintillation counting and standardized to 
cell protein. 

Statistical Analysis 

The analysis of the 1 H-NMR spectra was performed using an 
in-house written code in MATLAB® 7.1 software (MathWorks, 
Natick, MA) [38,44] and is described in detail in the Text SI. 
Briefly, spectra were segmented into bins and processed by 
principal component analysis (PGA) and by projection pursuit for 
exploratory data analysis (PPEDA) using Kurtosis as a projection 
index. Both types of analyses gave similar results regarding 
statistically significant differences between wildtype and NpcT 1 ' 
datasets. The results shown are from PPEDA. Following the initial 
exploratory data analysis, additional NMR spectra were acquired 
to achieve larger data sets for quantitative analysis of selected 
metabolites. 

The Shapiro-Wilk test was used to verify the normality of the 
distribution of all other data sets [45] . All NpcT 1 ' and age-matched 
wildtype control data sets satisfied the conditions of a normal data 
distribution according to the Shapiro-Wilk Test for thresholds 
p = 0.01, p = 0.05 and p = 0.1, with the exception of the qPCR 
analysis of mRNA levels of Gfap. A student's two-tailed t-test was 
used for the comparison of normally distributed data from NpcT 1 ' 
and age-matched wildtype control, and the Mann-Whitney U test 
was used for the comparison of mRNA levels of Gfap. Statistical 
significance was assumed as p<0.05. 

Results 

Metabolomic differences between wildtype and NpcV A 
brain regions in early stages of NPC disease 

To determine whether loss of NPC1 has consequences for 
energy metabolism in the brain, we first performed an unbiased 
metabolomics analysis of aqueous extracts of cerebellum, hippo- 
campus and cerebral cortex tissue of wildtype and NpcT / ' mice by 
1 H-NMR spectroscopy. Brain tissue was analyzed at 3, 4, 5, and 7 
weeks of age, representing pre-symptomatic (3 and 4 weeks), early 
symptomatic (5 weeks) and late stage (7 weeks) disease. Figure 1A 
shows a sample 1 H-NMR spectrum of an aqueous extract of 
wildtype cerebellum. Spectra for all samples were normalized to 
unit sum and mean-centred prior to multivariate exploratory 
analysis via projection pursuit (PPEDA). PPEDA is similar to the 
more commonly used principal component analysis (PGA), but is 
less sensitive to outliers than PCA [46-48]. PC A employs linear 
combinations of the original variables of high-dimensional data to 
obtain new, fewer variables (termed principal components; PC) 
that account for the most variance, i.e., the most information, in 
the original data. Typically, only a few PCs are sufficient for a 
good representation of the data. The projections of the data onto 
the new coordinate system defined by these PCs are depicted in so- 
called scores plots as shown in Figure IB. The scores plots allow an 
easy visual assessment of (dis)similarity among the data sets by the 
extent of clustering of data points derived from each data set. The 
scores plots resulting from our PCA analysis revealed significant 
differences between the metabolic profiles of wildtype and NpcT 1 ' 
cerebellum at 4, 5 and 7 weeks of age, but not at 3 weeks of age 
(Fig IB). In the cerebral cortex, scores plots showed significant 
differences only for 7-week old mice (Figure SI A). No significant 
differences between genotypes were found in the metabolomic 
analysis of hippocampal samples (Figure SIB). 
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Figure 1. Metabolomics analysis of wildtype and NpcT h cerebellum. 1 H-NMR spectroscopy of aqueous extracts from wildtype (WT) and 
Npcl cerebellum. A) Sample spectrum. Peaks: 1 ) Branched chain amino and ketoacids, 2) Lactate, 3) Alanine, 4) GABA, 5) Acetate/acetyl-CoA, 6) NAA, 
7) Glutamate, 8) Glutamine, 9) a-ketoglutarate, 10) Creatine, 1 1) Choline, 12) Taurine, and 13) myo-lnositol. B) PPEDA scores plots of sample sets of WT 
and Npcl cerebellum at 3-, 4-, 5- and 7- weeks of age. The separation of blue (WT) and red (Npcl~'~) data points indicates significant differences 
between WT and NpcV 1 ' sample sets at 4-, 5-, and 7- weeks of age. 
doi:1 0.1 371 /journal.pone.0082685.g001 



Increased lactate and decreased acetate/acetyl-CoA 
levels in the cerebellum of 3-week old NpcV A mice 

The scores plots alone can only show whether wildtype and 
Npcl''' brain regions differed in their overall metabolic profile, but 
they do not yield information about particular metabolites. On the 
other hand, the levels of specific metabolites could be different 
between the genotypes without causing a significant difference in 
the overall metabolic profile analyzed by PCA. We therefore 
deconvolved and integrated all spectra and determined the peak 
areas for all metabolites that were clearly distinguishable in the 
spectra. The peak areas were standardized to the total peak area to 
reflect the relative levels of energy metabolites in the aqueous 
extract. The most striking differences between wildtype and 
Npcl''' brain extracts were found in lactate and acetate/acetyl- 
CoA, which are the products of pyruvate reduction in anaerobic 
glycolysis and of oxidative decarboxylation of pyruvate during 
aerobic glycolysis, respectively. The resonance frequencies of the 



CH 3 -group of acetate and acetyl-CoA overlap near 1.92 ppm, and 
can therefore not be distinguished in the 'H-NMR spectra (Fig 1). 
Lactate levels were significantly higher in Npcl' 1 ' cerebellum and 
cerebral cortex compared to wildtype, and increased with disease 
progression in Npcl' 1 ' cerebellum (Fig 2 and Figure S2). In 
contrast, acetate/ acetyl-CoA levels were markedly decreased in 
the cerebellum and cerebral cortex of Npcl'*' mice aged 3 or 7 
weeks. In NpcT 1 ' hippocampus, acetate/acetyl-CoA levels were 
significandy decreased at both 3 and 5 weeks of age; however, 
lactate levels were unaltered (Figure S3). Alanine, which is 
catabolized to pyruvate by transamination, was increased in the 
cerebellum of 7-week old Npcf'~ mice (Fig 2). 

Our metabolomics analysis also revealed a decrease in the levels 
of the inhibitory neurotransmitter gamma-aminobutyric acid 
(GABA) in Npcl' 1 ' cerebellum at 5 and 7 weeks of age and in 
NpcT 1 ' cerebral cortex at 7 weeks of age, while glutamate levels 
were not altered (Fig 2 and Figure S2). Hippocampal Npcf ' 
samples showed a slight decrease in GABA levels at 3 weeks and 
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Figure 2. Altered energy metabolite levels in NpcT 1 ' cerebellum. Aqueous extracts of the cerebella from 3-, 4-, 5-, or 7-week old wildtype (WT) 
and NpcT 1 ' mice were analyzed by 'H-NIVIR spectroscopy. Spectra were deconvolved and integrated. Peak areas were standardized to total peak area. 
NAA: N-acetylaspartate. Inosit: myo-lnositol. Graphs in each column show data from the same set of WT and NpcV 1 ' mice of one age. The small 
numbers inside the bars for lactate indicate the number of mice in each group. Data are expressed as percent of the average of WT samples of the 
same age. Shown are the mean ± SEM. * p<0.05, NpcT 1 ' vs. WT. 
doi:10.1371/journal.pone.0082685.g002 



an increase in glutamate levels at 5 weeks of age (Figure S3). Thus, 
the neurotransmitter balance appeared to be shifted towards 
increased excitatory or decreased inhibitory neurotransmission in 
all brain regions analyzed. The resonance frequencies of glutamine 
and glutamate largely overlap, except for a region between 2.04 to 
2.07 ppm that corresponds mostly to glutamate. Peaks in this 
region were integrated for quantification of glutamate, and regions 
with overlapping resonance frequencies were quantified as the sum 
of glutamate and glutamine (Glu + Gin). There were no 
differences between genotypes in these peak ranges, suggesting 
that glutamine levels were not altered in NpcT 1 ' cerebellum (Fig 2). 

In spite of the known Purkinje cell loss in the cerebellum of 
older NpcT 1 ' mice, the levels of the neuron-derived metabolite N- 
acetylaspartate (NAA) were not different between genotypes, with 
only a non-significant (p = 0.08) decrease vaNpcT'' cerebellum at 7 
weeks of age (Fig 2). The levels of the osmolytes taurine and myo- 
inositol were slightly decreased in NpcT 1 ' cerebellum at 5 and 3 
weeks of age, respectively (Fig 2). The lack of alterations in myo- 
inositol content at more advanced stages of the disease was 
somewhat surprising, since the levels of the glia-derived myo- 
inositol commonly increase during glial activation and/or 
proliferation. Creatine levels remained unaltered in all samples 
(not shown). All brain areas of NpcT 1 ' mice 3 weeks of age and 
older showed a striking decrease in their levels of choline (Fig 2), 
which plays a key role in brain development as an important 
precursor of one-carbon metabolism and phospholipid synthesis 
[49]. 



Alterations in glucose and pyruvate metabolism in 
NpcV A cerebellum 

The increased lactate levels suggested that anaerobic glycolysis 
was increased in the NpcT 1 ' brains, while the decrease in acetate/ 
acetyl-CoA levels pointed to a decrease in the oxidative 
decarboxylation of pyruvate to acetyl-CoA for oxidative metab- 
olism of glucose. Given that there is limited fatty acid oxidation in 
the brain [50], the majority of acetyl-CoA found in the brain is 
derived from pyruvate. To investigate potential alterations in 
glucose metabolism in NPC 1 -deficient brain tissue, we first 
measured mRNA expression levels of several key enzymes of 
glycolysis in wildtype and NpcT ' cerebellum by qPCR. A 
schematic overview of glycolysis and its enzymes is shown in 
Fig 3A. The main rate-determining step of glycolysis is catalyzed 
by phosphofructokina.se- 1 (Pfk). At 3 weeks of age, we found no 
differences in the mRNA levels of two of the three PFK 1 isoforms 
[Pfkl and Pflcp), and decreased mRNA levels of the PJkm isoform in 
Npcl''' cerebellum compared to wildtype (Fig 3B). However, at 5 
weeks of age, mRNA levels of all three isoforms of PFK1 were 
markedly increased in NpcT 1 ' cerebellum (Fig 3B), suggesting that 
at this stage of the disease, glycolysis is upregulated on the gene 
expression level. In contrast, the mRNA levels of hexokinase 1 
(Hxkl) and pyruvate kinase type M (Pkni), which catalyze the two 
additional exergonic, highly regulated steps of glycolysis, were 
unchanged in NpcT 1 ' cerebellum at either 3 or 5 weeks of age 
(Fig 3C and D). Glucose-6-phosphate dehydrogenase (G6pd), 
which catalyzes the rate-limiting entrance of glucose-6-phosphate 
into the pentose-phosphate pathway, appeared decreased in 
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NpcT ' cerebellum at 3 weeks of age, but this decrease was not 
significant (p = 0.07; Fig 3E). To gain insight into the relative 
contributions of neurons and astrocytes, we measured mRNA 
levels of enolase, which is expressed in different isoforms in the two 
cell types [51]. While mRNA levels of astrocytic alpha-enolase 
{Enol) were increased in NpcT 1 ' cerebellum as early as 3 weeks of 
age, mRNA levels of neuron-specific enolase (Eno2) were unaltered 
at early stages of the disease and decreased in NpcT 1 ' cerebellum at 
7 weeks of age (Fig 3F and G). In addition, mRNA levels of the 
mostly astrocytic isoform of lactate dehydrogenase (Ldha), which 
catalyzes the production of lactate in the last step of anaerobic 
glycolysis, were significantly increased in NpcT 1 ' cerebellum at 5 
weeks of age (Fig 3H), in line with the increased lactate levels 
observed by NMR-spectroscopy (Fig 2). Together, these findings 
indicated differences in the transcriptional regulation of glucose 
metabolism, and a possible upregulation of glycolysis at the early 
symptomatic stage of the disease, in particular in astrocytes. 

Immunoblot analyses of HK 1 , PFK 1 and PKM protein levels in 
NpcT 1 ' cerebellum showed a pattern similar to the mRNA 
analysis. While HK1 levels were not altered (Fig 4A), protein 
levels of PFK1-P, which is highly expressed in cerebellum [52], 
were increased in NpcT 1 ' cerebellum at 5 weeks but not at 3 weeks 
of age (Fig 4B), further pointing to an increase in glycolysis in 
cerebellum at the early symptomatic stage of NPC disease. Total 
PKM protein levels were unaltered, however, the protein levels of 
the PKM2 isoform were increased in NpcT 1 ' cerebellum at 3 and 5 
weeks of age (Fig 4G). 

The observed increase in lactate levels, and increased levels of 
Enol and Ldha mRNA suggested increases in astrocytic glycolysis 
in NpcT 1 ' cerebellum. To investigate whether glycolysis may also 
be altered in NpcT 1 ' neurons, we measured the uptake of 
radiolabeled 2-deoxyglucose into primary wildtype and NpcT 1 ' 
cortical neurons after 9 days of culture. 2-Deoxyglucose is 
phosphorylated by hexokinase but not metabolized further, 
therefore its accumulation in the cells reflects glucose uptake and 
phosphorylation for glycolysis and pentose phosphate pathway 
[53]. NpcT 1 ' neurons had significantly higher 2-deoxyglucose 
uptake than wildtype primary neurons (Fig 4D), suggesting that 
glucose uptake and glycolysis can increase in NpcT 1 ' neurons. 

Decreased oxidative metabolism of pyruvate in NpcV A 
cerebellum 

Pyruvate derived from glycolysis can enter several different 
pathways, including the reversible transamination with alanine or 
the reduction to lactate. Pyruvate is also transported into 
mitochondria, where it can undergo carboxylation to oxaloacetate 
or oxidative decarboxylation to acetyl-CoA, which can itself be 
further oxidized in the citric acid cycle (schematic overview shown 
in Fig 5A). The increased lactate to acetate/acetyl-CoA ratio 
observed by NMR spectroscopy in NpcT 1 ' brain at all ages 
suggested that mitochondrial oxidative metabolism of pyruvate 
was decreased, and that the cytosolic reduction of pyruvate to 
lactate was increased to supply NAD + for glycolysis. 

To investigate potential alterations in the mitochondrial 
metabolism of pyruvate in NpcT 1 ' cerebellum, we measured the 
mRNA levels of several key enzymes of pyruvate metabolism and 
of the malate aspartate transporter Aralar 1 , which is part of a 
mitochondrial shuttle for the indirect transfer of NAD into the 
cytosol. The oxidative decarboxylation of pyruvate to acetyl-CoA 
is catalyzed by the pyruvate dehydrogenase complex PDH in the 
mitochondrial matrix. Expression analysis of PDH showed lower 
mRNA levels of the Pdha subunit in NpcT 1 ' cerebellum at 3 weeks 
of age (Fig 5B). Moreover, immunoblot analysis revealed a 
decrease in the total levels of PDH in NpcT 1 ' cerebellum at 5 weeks 



of age, and an increase in the levels of phosphorylated PDH at 3 
and 5 weeks of age (Fig 5C). PDH is inhibited by phosphorylation; 
therefore, these findings strongly suggested a decrease in the 
oxidation of pyruvate to acetyl-CoA, in line with the decreased 
acetate/acetyl-CoA levels observed by NMR spectroscopy (Fig 2). 
The increased mRNA levels of PDH kinases 1 and 2 in NpcT 1 ' 
cerebellum at 5 weeks of age (Fig 5D and E; Pdkl and Pdk2) 
suggested a longer-term increase in PDH phosphorylation on the 
gene expression level in Npcl~'~ cerebellum at this age. Given that 
PDK1 and PDK2 are more prevalent in neurons and astrocytes, 
respectively [54], these alterations in PDH likely occurred in both 
cell types. Cerebellar mRNA levels of pyruvate carboxylase (Pcx), 
which converts pyruvate to oxaloacetate only in astrocytes, were 
increased in 3- and 5-week old NpcT 1 ' mice (Fig 5F), while mRNA 
levels of cytosolic malic enzyme (Mel), which catalyzes the 
interconversion of pyruvate and malate mainly in the direction 
of pyruvate synthesis, were decreased in NpcT 1 ' cerebellum at 3 
weeks of age (Fig 5G). Mitochondrial malic enzyme 2 (Me2) 
mRNA levels were also decreased in NpcT 1 ' cerebellum at 3 weeks 
of age (Fig 5H), while mRNA levels of the mainly neuronal, 
mitochondrial malic enzyme 3 (Me3) and the malate aspartate 
shuttle aralar (AralarT) were unchanged, except at 7 weeks of age 
(Fig 51 and J). The decreased mRNA levels of neuronal Eno2 and 
Aralarl in NpcT 1 ' cerebellum at 7 weeks of age (Fig 3G and 5J) 
likely reflected Purkinje neuron death, as mRNA levels of the 
Purkinje cell marker calbindin (Calbl) and of the neuron-specific 
neurofilament heavy peptide (Nejh) were decreased by 60% and 
25% respectively at this age (Fig 5K and L). 

Together, these findings demonstrated that alterations in 
cerebellar glucose metabolism (upregulation of Enol and Pcx 
mRNA, higher PKM2 protein, decreased Pdha, Mel and Me2 
mRNA levels, increased lactate and decreased acetate/acetyl-CoA 
levels) developed early in NPC disease, and could already be 
detected in 3-week old, pre-symptomatic NpcT 1 ' mice. Dysregu- 
lation of gene expression and of the metabolic profile of Npcl ' 
cerebellum was even more pronounced at 5 weeks, when 
symptoms become apparent in NpcT 1 ' mice. 

Mitochondrial alterations in Npcl ~'~ brain 

Next, we asked whether the upregulation of glycolytic enzymes 
and the increased lactate production in NpcT 1 ' cerebellum were a 
response to decreased mitochondrial abundance and/or mito- 
chondrial defects, either of which would lower glucose oxidation 
and decrease mitochondrial ATP production. A key regulator of 
mitochondrial biogenesis and function is the peroxisome pro- 
liferator-activated receptor gamma coactivator 1 alpha (PGCTot), 
which interacts with a variety of transcription factors involved in 
energy metabolism and mitochondrial function to increase 
mitochondrial biogenesis. PGC-la mRNA (Ppargcld) levels were 
not altered in NpcT*' cerebellum at 3 weeks, but were markedly 
increased compared to wildtype at 5 weeks of age, and decreased 
at 7 weeks of age (Fig 6A), suggesting that at 5 weeks of age, the 
cellular response to increase mitochondrial biogenesis was 
activated in NpcT 1 ' cerebellum. We therefore determined the 
levels of a mitochondrial marker protein, the voltage dependent 
anion channel (VDAC1), and levels of mitochondrial DNA as 
measures for mitochondrial abundance. While VDAC1 protein 
levels remained unaltered throughout disease progression, mito- 
chondrial DNA levels were substantially decreased in the cerebella 
of 5-week old Npcl' 1 ' mice (Fig 6B and C), suggesting a loss of 
normal mitochondria even though Ppargcla mRNA levels were 
increased. Additional alterations in mitochondria were seen in a 
decrease in mRNA levels of the mitochondria-encoded subunit 2 
of cytochrome c oxidase (Mtco2) and a marked increase of mRNA 
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Figure 3. Increased expression of glycolytic enzymes in NpcT 1 ' cerebellum at 5 weeks of age. A) Schematic representation of glycolysis 
and connected pathways showing mRNA (italics) and proteins analyzed in this study. PPP: pentose phosphate pathway. B - H) Cerebellar RNA 
extracts were prepared from 3-, 5- or 7-week old wildtype (WT) and NpcV 1 ' mice, and mRNA levels of the target genes were analyzed by qPCR using 
cyclophilin (Ppia) as housekeeping gene. B) Phosphofructokinase isoforms (Pfk!, Pfkm, Pfkp), Q Hexokinase 1 (Hxk1); D) Pyruvate Kinase Type M (Pkm), 
E) astrocytic alpha-Enolase (Eno1), F) neuronal Enolase (Eno2), G) Lactate Dehydrogenase (Ldha), H) Glucose-6-phosphate dehydrogenase (G6pd). Data 
are shown as mean ± SEM. * p<0.05, NpcV'~ vs. WT. 
doi:1 0.1 371 /journal.pone.0082685.g003 



levels of uncoupling protein 2 (Ucp2) in NpcT ' cerebellum at all 
ages (Fig 6D and E). These observations suggested alterations in 
mitochondrial abundance, quality and electron transport chain 
components that could affect energy homeostasis in NpcT 1 ' 
cerebellum. 

To investigate whether oxidative phosphorylation may be 
impaired in NpcT 1 ' brain, we measured ATP production by 
mitochondria isolated from wildtype and NpcT f ~ brain at 3 and 7 
weeks of age. Mitochondrial polarization was verified with 
acridine orange staining (Fig 6F). In the presence of pyruvate 
and malate as energy substrates, mitochondria isolated from 
NpcT f ~ brain at 3 weeks produced less ATP than mitochondria 
from wildtype brain, in line with a decrease in PDH activity and 



impaired oxidative phosphorylation in NpcT ' brain. On the other 
hand, in the presence of succinate, ATP generation was 
unchanged (Fig 6G), which may suggest that most components 
of the electron transport chain remained relatively unaffected in 
mitochondria from NpcT 1 ' brain. Interestingly, mitochondria 
isolated from wildtype and NpcT 1 ' brain at 7 weeks of age showed 
no differences in ATP generation (Fig 6H). While we cannot 
explain why isolated mitochondria from NpcT 1 ' brain at later 
stages of the disease appear to be normal in this assay, the findings 
reflect the lack of differences between genotypes at 7 weeks of age 
in most of our metabolic gene expression analyses (Fig 3 and 5). It 
is important to note, however, that assays using isolated 
mitochondria cannot fully reflect mitochondrial function in the 
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Figure 4. Increased phosphofructokinase levels in Npcr 1 ' cerebellum at the early symptomatic stage of the disease. A to C) Cerebellar 
tissue homogenates were prepared from 3- or 5-week old wildtype (WT) and Npc1~'~ mice and analyzed by immunoblotting using antibodies against 
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brain. Isolation procedures and the assay environment strongly 
influence the measurements, and alterations in mitochondrial 
function in the brain are not necessarily recapitulated by isolated 
mitochondria. 

To further investigate energy homeostasis in NpcT 1 ' cerebellum, 
we determined the levels of phosphorylated AMP-activated kinase 



(AMPK) by immunoblotting. A drastic loss of oxidative phos- 
phorylation capacity would be expected to have a significant 
impact on energy homeostasis and likely lead to phosphorylation 
and activation of AMPK. However, levels of phosphorylated 
AMPK were unchanged in NpcT 1 ' cerebellum at 3, 4, and 5 weeks 
of age (Fig 61), suggesting that there was no significant overall 
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Figure 5. Altered expression of pyruvate metabolism-related genes in Npcr 1 ' cerebellum. Cerebellar tissue homogenates and RNA were 
prepared from 3-, 5-, or 7-week old wildtype (WT) and NpcV'~ mice. A) Schematic representation of pathways of pyruvate metabolism showing mRNA 
(italics) and proteins analyzed in this study. The grey shading represents the mitochondrial compartment. B) mRNA levels of pyruvate dehydrogenase 
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Figure 6. Mitochondrial alterations in NpcT 1 ' cerebellum. A) Cerebellar RNA was prepared from 3-, 5- or 7-week old WT and Npc1~'~ mice. 
Target gene mRNA levels were measured by qPCR using primers against PGC1 a {Ppargda) standardized to Ppia. B) Immunoblots of cerebellar tissue 
homogenates from 3- or 5-week old wildtype (WT) and Npc1~'~ mice probed with antibodies against VDAC1 and actin as loading control. C) Genomic 
DNA was prepared from cerebella of 5-week old WT and Npcl mice. Mitochondrial DNA was determined by qPCR using primers against 
mitochondrial cytochrome c oxidase subunit 1 (Mtcol) DNA per nuclear Ndufvl DNA, and standardized to WT. D and E) Cerebellar RNA was prepared 
as in panel A. Target gene mRNA levels were measured by qPCR using primers against mitochondrially encoded subunit 2 of cytochrome c oxidase 
(Mfco2) or uncoupling protein 2 {Ucp2) standardized to Ppia. F to H) Mitochondria were isolated by differential centrifugation from combined 
homogenates of the cerebella, hippocampi and cerebral cortices of wildtype (WT) and Npc1~'~ mice at 3 and 7 weeks of age. F) Mitochondria were 
stained with nonyl acridine orange to verify the maintenance of their membrane potential. G and H) The rate of ATP generation was measured in the 
presence of 1 mM pyruvate, 1 mM malate and 300 uM ADP (Pyr/Mal) or in the presence of 5 mM succinate and 300 uM ADP as luminescence 
generated by the ATP-dependent D-luciferin/luciferase reaction as luminescence per second per protein. F) Mitochondria isolated from 3-week old 
mice. G) Mitochondria isolated from 7-week old mice. Data shown in panels F - H) are expressed as percent of WT of the same experiment, and are 
derived from three independent experiments each with two independent mitochondria preparations of each genotype and each age. I) Cerebellar 
tissue homogenates were prepared from 3-, 4- or 5-week old WT and NpcV'~ mice and analyzed by immunoblotting using antibodies against 
phospho-AMPKa (Thr172), AMPKoc, phospho-AMPKB (Ser108), AMPKB, and actin as loading control. Bar graphs show the ratio of band intensities of 
phosphorylated to total AMPK expressed as percent of the average of WT on the same membrane. All bar graphs show the mean ± SEM. * p<0.05, 
Npcf'- vs. WT. 
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energy deficiency, possibly due to increased glycolytic ATP 
production. 

Activation of antioxidant defense systems in NpcV A 
cerebellum 

Next, we measured mRNA levels of several transcription factors 
or known target genes associated with stress responses. Glial 
fibrillary acidic protein [Gfap), a marker for astrocyte activation, 
was significandy increased in Npcl' 1 ' brain as early as 3 weeks of 
age (Fig 7A), in accordance with previous reports [55-57]. 
Glycolytic gene expression and anaerobic glycolysis are also 
induced by activation of hypoxia-induced factor la (HIF-la, 
Hi/la) [58]. However, Hifla mRNA levels were not different 
between genotypes (Fig 7B), and we were unable to reliably detect 
HIF-la by immunoblotting in either wildtype or Npcl' 1 ' brain 
homogenates. In contrast, mRNA levels of the nuclear factor E2 
related factor 2 (Nrf2, Nfe2l2), which is a key element of a 
protective response to oxidative stress, and of one of its target 
genes, heme oxygenase 1 (HOI, Hmoxl), were markedly higher in 
Npcl''' cerebellum beginning at 3 weeks of age, and increased with 
disease progression (Fig 7C and D). The early activation of 
antioxidant responses was also apparent in the increased mRNA 
levels of superoxide dismutase 1 (Sodl) and the increased protein 
levels of superoxide dismutase 2 (MnSOD) in Npcl' 1 ' compared to 
wildtype cerebellum (Fig 7E and F). Moreover, we observed higher 
levels of phosphorylation of both isoforms of the c-jun N-terminal 
kinase (JNK) in NpcT 1 ' cerebellum at 3 and 5 weeks of age (Fig 7G). 
Phosphorylation and activation of JNK is a common response to 
oxidative stress and a wide range of other stress signals, suggesting 
that oxidative stress increases very early during NPC disease 
pathogenesis. 

Discussion 

Here we have investigated brain energy metabolism in a mouse 
model of NPC disease. Using an unbiased metabolomics strategy 
and targeted gene expression analyses of wildtype and NpcT 1 ' 
brain tissue, we found i) significant alterations in glucose and 
pyruvate metabolism that suggested a pre-symptomatic decrease in 
the oxidative metabolism of glucose, and a compensatory increase 
in glycolysis with disease progression, and ii) a pre-symptomatic 
activation of astrocytic antioxidant response systems, which 
suggests that oxidative stress plays a key role early in NPC disease 
pathogenesis. A graphic overview of the gene expression changes 
found in NpcT 1 ' cerebellum at 3 and 5 weeks of age is shown in 
Figure 8. 



Glycolysis and pyruvate oxidation in Npcl cerebellum 

The most striking alterations found in our metabolomics 
analysis were an increase in lactate and a decrease in acetate/ 
acetyl-CoA levels in NpcT 1 ' cerebellum and cerebral cortex as 
early as 3 weeks of age (Fig 2). The decrease in acetate/acetyl-CoA 
pointed to an impaired oxidative decarboxylation of pyruvate to 
acetyl-CoA for further oxidation in the citric acid cycle. High 
lactate levels are often indicative of increased anaerobic glycolysis 
due to hypoxia or mitochondrial dysfunction, and have been 
observed in brain and cerebrospinal fluid in cerebral ischemia, 
during physiological aging, and in certain models of neurodegen- 
erative diseases associated with mitochondrial dysfunction, such as 
Huntington disease [59-63]. In addition, lactate levels can 
increase transiently with increased neuronal activity [64—66] . 

The increased mRNA levels oiLdha, Enol, and PJk isoforms and 
the higher PFK1 protein levels observed in Npcl''' cerebellum at 5 
weeks of age indicated that glycolysis was upregulated and 
contributed to the higher lactate levels in NpcT 1 ' brain at this 
age (Fig 3 and 4). Similarly, published microarray data from 
NpcT 1 ' mouse liver show increased expression of several glycolytic 
enzymes, including PJk and Ldh isoforms [67]. Given that Ldha and 
Enol are predominantly expressed in glia, which are also the main 
glycolytic cell type and the main source of lactate in normal brain, 
the increased expression of these isoforms suggested an increase in 
glycolysis in Npcl''' astrocytes. The increased 2-deoxyglucose 
uptake into primary Npcl''' neurons (Fig 4D) suggested that 
glycolysis can also increase in neurons in Npcl''' brain. Whereas in 
astrocytes, increased glycolysis may serve to provide lactate and 
other energy substrates to neurons [31,68,69], an upregulation of 
neuronal glycolysis can have detrimental effects by limiting the 
availability of glucose for the pentose phosphate pathway and thus 
for synthesis of NADPH and regeneration of glutathione [70]. 

Although lactate levels were increased at all stages of NPC 
disease (Fig 2), the upregulation of glycolysis in Npcl''' cerebellum, 
at least on gene and protein expression level, appeared to occur 
during progression to the early symptomatic stage rather than pre- 
symptomatically. Thus, at 3 weeks of age, only Enol mRNA levels 
were increased while Pfkm mRNA levels were even slighdy 
decreased in Npcl''' cerebellum (Fig 3). With progression from 
early symptomatic to late stage NPC disease, glycolytic gene 
expression in NpcT 1 ' cerebellum appeared to normalize, even 
though lactate levels remained elevated (Fig 3). A recent 
proteomics study using the same NpcT 1 ' mice found decreased 
levels of glycolytic PKM1/2 and PFK-M in Npcl' 1 ' cerebellum at 
1 and 3 weeks of age respectively [71], in agreement with our 
findings at 3 weeks. However, whereas our results suggest that 
glycolysis increases in the early stages of NPC disease, these 
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hypoxia-inducible factor 1 alpha (Hifla), nuclear factor erythroid 2-related factor 2 (Nfe2l2), heme oxygenase 1 (Hmoxl), and superoxide dismutase 1 
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authors concluded that glycolysis was generally downregulated in 
NPC disease, based on their findings of decreased levels of 
glycolytic triosephosphate isomerase in NpcT 1 ' cerebellum at 5 
weeks of age, and of decreased deoxyglucose uptake in cultured 
primary fibroblasts from human NPC patients [71]. In vivo 
measurements of brain glucose uptake by 2-Deoxy-2-[ 1!! F]fluoro- 
d-glucose positron emission tomography have shown glucose 
hypometabolism in one patient with adult onset NPC disease [72], 
whereas two other case reports of NPC patients examined with the 
same technique found unaltered or increased glucose metabolism 
in the cerebellum and hypometabolism only in the cerebral cortex 
[73,74]. Clearly, NPC 1 -deficiency alters glucose metabolism in a 
complex manner that is highly dependent on the stage of the 
disease and on the brain region examined. Further investigations 
by in vivo measurements of metabolic flux or by neuroimaging 
approaches are needed to clarify the role of glycolytic activity in 
the NPC 1 -deficient brain [75,76]. 

We also observed an increase in the levels of the M2 isoform of 
pyruvate kinase (PKM2, Fig 4C) in Npd''' cerebellum. PKM2 is 
commonly expressed in embryonic and cancer tissue, and, in 



contrast to PKM1, it is under complex allosteric regulation to 
control the flux of glycolytic intermediates into biosynthetic 
pathways [77,78]. The consequences of increased PKJVI2 levels 
in NpcT 1 ' cerebellum are unclear, but expression of an embryonic 
isozyme may point to disturbances in cell cycle or differentiation. 
Similarly, higher mRNA levels of the embryonic aldolase B have 
been found in Npd''' murine liver [67]. 

The increased lactate and decreased acetate/ acetyl-CoA levels 
in NpcT 1 ' cerebellum at 3 weeks of age without a clear 
upregulation of glycolysis at this age (Fig 2 and 3) suggested that 
pyruvate metabolism was affected by NPC 1 -deficiency. Alterations 
in pyruvate metabolism were apparent in the markedly increased 
Pcx and decreased Pdha, Mel and Me2 mRNA levels in Npd''' 
cerebellum at 3 weeks of age (Fig 5). The strong trend to increased 
inhibitory phosphorylation of PDH in NpcT 1 ' cerebellum at 3 
weeks of age, the significantly increased phosphorylated PDH and 
the decreased total PDH protein in NpcT 1 ' cerebellum at 5 weeks 
of age (Fig 5) further indicated that an impaired oxidative 
decarboxylation of pyruvate to acetyl-CoA contributed to the high 
lactate and low acetate/acetyl-CoA levels in NpcT 1 ' cerebellum. 
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Figure 8. Overview of the alterations in gene and protein expression in NpcT 1 ' cerebellum and the proposed effects on metabolic 
pathways. Schematics show glycolysis and pyruvate metabolism annotated with the mRNAs and proteins analyzed in this study in cerebellum from 
3- or 5-week old WT and Npd~'~ mice. Green squares: increased mRNA or protein levels. Red ovals: decreased mRNA or protein levels. No colour: 
unchanged mRNA or protein levels. PPP: pentose phosphate pathway. The grey shading represents the mitochondrial compartment. The proposed 
effects on the metabolic pathways are depicted by the thickness of the arrows. 
doi:1 0.1 371 /journal. pone.0082685.g008 
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Impaired PDH activity could also contribute to the decreased ATP 
generation in the presence of pyruvate/malate that we observed in 
assays with mitochondria isolated from NpcT 1 ' brain at 3 weeks of 
age (Fig 6). In a previous study, Cologna et al. also found 
significandy decreased PDH protein levels and a strong trend to 
increased pyruvate levels in the cerebella of 1-week old NpcT 1 ' 
mice [71]. Other reports have shown decreased mRNA expression 
of PDH isoforms in NpcT 1 ' cerebellum at a late stage of the 
disease, as well as in NpcT 1 ' liver [27,67]. Moreover, increased 
mRNA levels of different PDH Kinase isoforms have been 
reported for liver samples from Npcl 1 mice aged 1 week and 
older, and for NPC 1 -deficient fibroblasts [67,79]. Thus, a decline 
in PDH activity seems to be an early and persistent consequence of 
NPC 1 -deficiency in many cell types. Our finding that both Pdkl 
and Pdk2 mRNA levels were increased in NpcT 1 ' cerebellum (Fig 5) 
suggested that PDH phosphorylation increased in both neurons 
and astrocytes, since neurons express PDK1, and PDK2 is 
predominantly present in astrocytes [54]. The upregulation of 
astrocytic Pcx in NpcT 1 ' cerebellum at 3 weeks of age and 
concomitant downregulation of Mel and Me2 (Fig 5) may serve to 
counteract a build-up of pyruvate and the subsequent feedback 
inhibition of glycolysis, as has been described in other cell types 
[80]. 

The brain is highly sensitive to a decrease in PDH activity 
because of its strong reliance on the oxidative metabolism of 
glucose [28]. Congenital PDH deficiencies manifest with neuro- 
logical symptoms such as ataxia and seizures [81], and decreases in 
PDH activity are also observed in Alzheimer disease and 
physiological aging [82,83]. In astrocytes, PDH activity is 
normally kept low by high basal phosphorylation [54], and 
pyruvate is diverted into lactate formation and anaplerotic 
pyruvate carboxylation [84]. In neurons, however, maintenance 
of dephosphorylated, active PDH is essential to achieve the most 
energy-efficient complete oxidation of glucose [54]. Impaired 
oxidative decarboxylation of pyruvate also prevents the use of 
lactate or alanine as energy substrates, leaving increased glycolysis 
or oxidation of glutamine as the main alternative energy sources. 
Thus, a decrease in PDH activity could lead to impaired energy 
generation and/ or to metabolic adaptations with possible conse- 
quences for neuronal redox status and neurotransmitter metabo- 
lism [85] . In addition, differences in neurotransmitter metabolism 
can also contribute to alterations in energy metabolism. Increased 
neuronal activity increases the neuronal energy requirement to 
maintain the membrane potential, and is commonly coupled to 
increased astrocytic lactate production [64,65]. If the imbalance in 
neurotransmitter levels and in particular the decrease in inhibitory 
GABA that was observed in Npcl''' cerebellum at 5 and 7 weeks of 
age (Fig 2) led to changes in neuronal activity, it could thus 
contribute to the increase in lactate production at later stages of 
the disease. On the other hand, in Npcl''' hippocampus, GABA 
levels were decreased at 3 weeks of age without a concomitant 
increase in lactate (Figure S3). Moreover, a variety of homeostatic 
synaptic scaling mechanisms can stabilize overall neuronal firing 
through compensatory downregulation of excitatory activity [86- 
88]. 

Other factors could also contribute to impaired mitochondrial 
energy generation in NpcT 1 ' brain. Our findings of decreased 
Mtco2 mRNA and increased Ucp2 mRNA levels in NpcT 1 ' 
cerebellum at all ages (Fig 6) could indicate defects in cytochrome 
c oxidase /complex IV, and possibly increased uncoupling of the 
electron transport chain, even though the role of UCP2 in the 
brain is not quite clear. Cologna et al. found decreased ATP 
synthase levels in Npcl' 1 ' cerebellum at 1 week of age, but no 
changes in ATP synthase levels at 3 weeks of age [71]. Using 



isolated brain mitochondria from 3-week old Npcf ' mice, we 
found a decrease in ATP production with pyruvate and malate as 
energy substrates, but no differences in the presence of succinate. 
Even with the caveat that assays in isolated mitochondria are 
limited by their artificial environment, these findings suggest that 
impaired oxidative metabolism of pyruvate may be more 
prominent than electron transport chain deficiencies. 

In several other neurodegenerative diseases associated with 
mitochondrial dysfunction, deficiencies in mitochondrial ATP 
production are associated with increased phosphorylation and 
activation of AMPK [89—93]. However, we observed no 
significant alterations in the levels of phosphorylated AMPK in 
cerebellum from 3-, 4- or 5-week old Npcl''' mice (Fig 6), 
suggesting that any deficiencies in oxidative phosphorylation in 
Npcl''' cerebellum were mild enough to be offset by compensatory 
metabolic adaptations. Such compensatory mechanisms could 
include the upregulation of Pcx in astrocytes in Npcl'*' cerebellum 
at 3 weeks of age and the upregulation of glycolytic gene 
expression (PFK, PJk, Ldhd) at 5 weeks of age (Fig 3 - 5), which 
could serve to counteract pyruvate build-up and energy deficits 
during the early symptomatic stage of NPC disease. Moreover, a 
proteomics analysis showed increased ATP synthase levels in 
NpcT 1 ' cerebellum at 5 weeks of age [71]. Our finding of increased 
mRNA levels of PGC-la, a transcriptional coactivator for 
mitochondrial biogenesis, in NpcT 1 ' cerebellum at 5 weeks of 
age (Fig 6) suggested a compensatory upregulation of mitochon- 
drial biogenesis, however, the lower levels of mitochondrial DNA 
observed in NpcT 1 ' cerebellum at the same age (Fig 6) indicated a 
decreased mitochondrial number or mitochondrial defects. A 
possible explanation for this seeming contradiction may come 
from recent work by Ordonez et al. showing an accumulation of 
depolarized, dysfunctional mitochondria due to defective mito- 
phagy in NPC 1 -deficient human embryonic stem-cell derived 
neurons [24]. We speculate that with progression to late stage 
NPC disease, compensatory mechanisms begin to fail. Thus, at 7 
weeks of age, PGC-la mRNA levels in NpcT 1 ' cerebellum were 
decreased (Fig 6) and metabolic gene expression appeared mosdy 
unaltered or decreased due to neuronal loss (Fig 3 and 5; Eno2, 
Aralarl). The lack of significant defects in ATP generation by 
mitochondria isolated from the brains of 7-week old Npcf'' mice 
(this study, Fig 6), but decreased ATP generation by mitochondria 
isolated from brain homogenates of 9-week old Npcf'' mice 
observed by Yu et al. [22], also indicate a progressive deterioration 
of mitochondrial quality. 

Early increase in oxidative stress and activation of 
antioxidant response systems 

The most striking pre-symptomatic alterations in NpcT 1 ' 
cerebellum in mRNA and protein expression were related to 
oxidative stress and antioxidant response systems, for example 
increases in the levels of Nfe2L2, Hmoxl, Sodl, and Ucp2 mRNA 
and of pJNK and MnSOD protein (Fig 7). Increases in markers for 
oxidative stress have been found previously in several cell and 
animal models of NPC, as well as in human NPC patients 
[8,25,27,67,94,95]. However, whereas our findings strongly 
suggest a role for oxidative stress early during NPC disease 
pathogenesis, a previous study by Cluzeau et al. had concluded 
that oxidative stress was a late occurrence in NPC disease, based 
on microarray data of liver taken from NpcT*' mice between the 
ages of 1 and 1 1 weeks [67]. In that study, pathway maps related 
to oxidative stress showed significant alterations in NpcT 1 ' liver 
only at 1 1 weeks of age [67]. Even if the difference to our results 
may in part be due to methodological differences, the much later 
onset of oxidative stress in NpcT 1 ' liver compared to brain seems to 
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suggest that in NPC disease, the cerebellum may be exposed to 
greater oxidative stress or be more sensitive to oxidative stress than 
the liver. Differences in the susceptibility to oxidative stress 
between mitochondria from liver and brain have been reported 
previously [96,97]. A recent study investigated the effects of 
supplementing Npcl''' mice with the antioxidant N-acetylcysteine 
beginning at 4 or 6 weeks of age, and noted a modest reduction in 
oxidative stress markers and a slight improvement in certain 
aspects of the NPC disease phenotype in the N-acetylcysteine- 
treated animals [98] . Our findings indicate that oxidative stress in 
the Npcl ' cerebellum increases before the onset of overt 
neurological symptoms, and suggest that early antioxidant 
treatment may be beneficial in NPC disease. 

Oxidative stress and mitochondrial dysfunction form a vicious 
cycle observed in many neurodegenerative disorders, where 
alterations in the electron transport chain and in oxidative 
phosphorylation lead to increased production of reactive oxygen 
species (ROS), which in turn promote further mitochondrial and 
metabolic dysfunction through oxidative modification and accel- 
erated degradation of glycolytic and mitochondrial proteins [99] . 
Interestingly, PDH appears to be highly sensitive to oxidative 
modification, as seen following cerebral ischemia or in acutely 
respiring yeast [99-101]. JNK activation by oxidative stress can 
also decrease PDH activity. Thus, in cortical neurons, activated 
JNK translocates to the mitochondrial outer membrane and 
initiates a phosphorylation cascade that results in inhibitory PDH 
phosphorylation [83,102]. In addition, JNK activation can 
increase PDK1 transcription [103]. 

The initial cause of a deterioration in mitochondrial quality in 
Npcl''' brain is unclear. One possibility is that impaired mitophagy 
leads to the accumulation of dysfunctional mitochondria that 
produce more ROS [104]. Several groups, including ours, have 
also shown that cholesterol trafficking defects in NPC 1 -deficient 
cells can lead to an increase of cholesterol in mitochondria, which 
in turn may affect their function [20-23]. An upregulation of 
glycolysis in Npcl''' neurons and/or decreased ME1 activity may 
further increase oxidative stress by decreasing NADPH synthesis 
for glutathione regeneration [70]. Initial increases in oxidative 
stress or initial metabolic alterations can thus cause further 
changes in energy metabolism and increase oxidative damage in 
NpcT 1 ' brain, leading to a cascade of subsequent adaptations as 
seen in many neurodegenerative diseases. 

Our study also provides further evidence that both neurons and 
astrocytes are affected by NPC 1 -deficiency, which is in agreement 
with previous reports of early microglial and astrocyte activation 
[13,57,105,106]. However, the observed pre-symptomatic meta- 
bolic alterations in astrocytes also raise questions as to the role of 
astrocytes in NPC pathology. On one hand, studies using chimeric 
mice [107], mice with cell-type specific deletion of NPC1 [108] or 
with cell-type specific expression of NPC1 in an NPC 1 -null 
background [109] have clearly demonstrated a neuron-inherent, 
glia-independent mechanism of Purkinje cell death in NPC1- 
deficient cerebellum. Moreover, depletion or expression of NPC 1 
in astrocytes alone neither caused nor rescued the NPC 
neurological phenotype [108,109], with the exception of one 
report [110]. On the other hand, given the close interdependence 
of neurons and astrocytes, it seems unlikely that astrocytes could 
show such striking and early metabolic alterations, and yet have no 
impact on NPC disease progression. Our findings cannot resolve 
this quandary; however, they clearly show that many of the 
alterations in NpcT 1 ' astrocytes can be neuroprotective, including 
the activation of Nrf2 and upregulation of Hmoxl (HOI) [Hi- 
ll 4]. Nrf2 activity increases astrocytic synthesis and secretion of 
glutathione [111,112], which, after hydrolysis, generates cysteine 



and glycine for neuronal glutathione synthesis [111,112]. HO-1 
also plays a mostly cytoprotective role, but it can also exacerbate 
neuronal and astrocytic injury after prolonged increased expres- 
sion [115,116], possibly contributing to the further deterioration 
seen at late stage NPC disease. For the development of additional 
therapeutic intervention strategies aimed at energy metabolism, it 
will be important to determine in more detail where and when 
different metabolic alterations develop, and whether these changes 
are potentially neuroprotective or part of the neuronal decline. 
Moreover, the significant alterations we observed prior to onset of 
symptoms once again underline the importance of early diagnosis 
and intervention. 

Supporting Information 

Figure SI PPEDA scores plots for the 'H-NMR spectra 
of wildtype and Npcl''' cerebral cortex and hippocam- 
pus. Scores plots representing the comparisons of wildtype (WT, 
blue squares) and NpcT*' (red squares) samples at 3, 4, 5 and over 
7 weeks of age as indicated. A) Cerebral Cortex samples. B) 
Hippocampal samples. Only sample sets of WT and Npcl''' 
cerebral cortices taken at 7 weeks of age were statistically 
significandy different. 
(TIF) 

Figure S2 Alterations in energy metabolite levels in 
Npcl''' cerebral cortex. Aqueous extracts of the cerebral 
cortices from 3-, 4-, 5-, and 7-week old wildtype (WT) and NpcT 1 ' 
mice were analyzed by 1 H-NMR spectroscopy. Spectra were 
deconvolved and integrated. Peak areas were standardized to total 
peak area. NAA: N-acetylaspartate. Inosit: myo-Inositol. Graphs 
in each column show data from the same set of WT and NpcT 1 ' 
mice of one age. The small numbers inside the bars for lactate 
indicate the number of mice in each group. Data are expressed as 
percent of the average of WT samples of the same age. Shown are 
the mean ± SEM. * p<0.05, Npcl''' vs. WT. 
(TIF) 

Figure S3 Alterations in energy metabolite levels in 
Npcl''' hippocampus. Aqueous extracts of the hippocampi 
from 3- and 5-week old wildtype (WT) and Npcl''' mice were 
analyzed by 1 H-NMR spectroscopy. Spectra were deconvolved 
and integrated. Peak areas were standardized to total peak area. 
NAA: N-acetylaspartate. Inosit: myo-Inositol. Graphs in each 
column show data from the same set of WT and Npcl''' mice of 
one age. The small numbers inside the bars for lactate indicate the 
number of mice in each group. Data are expressed as percent of 
the average of WT samples of the same age. Shown are the mean 
± SEM. * p<0.05, NpcT'' vs. WT. 
(TIF) 

Text SI Supporting methods. Detailed description of the 
acquisition of 1 H-NMR spectra and the statistical analysis by 
PPEDA. 
(PDF) 

Table SI Primer pairs used for qPCR analysis. Listed are 
the sequences of the primer pairs and the annealing temperatures 
used for qPCR analysis as well as the accession numbers of the 
genes of interest. 
(PDF) 

Table S2 Relative mRNA expression of potential house- 
keeping genes Ppia, Actin, and Rpll3a. The mRNA levels of 
each of the three potential housekeeping genes Ppia, Actin, and 
Rpll3a were measured in wildtype (WT) and NpcT 1 ' cerebellum 
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and their relative expression was calculated by the Pfaffl method 

with each of the other two as housekeeping genes. 

(PDF) 

Table S3 Primary antibodies used in this study. Listed 
are the sources for the primary antibodies used for immunoblot- 
ting. 
(PDF) 
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